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ABSTRACT 

We have conducted a survey of the CCS Jjv — 32 — 2i line toward 11 dark 
clouds and star-forming regions at 30" spatial resolution and 0.054 kms -1 
velocity resolution. CCS was only detected in quiescent clouds, not in active 
star-forming regions. The CCS distribution shows remarkable clumpy structure, 
and 25 clumps are identified in 7 clouds. Seven clumps with extremely narrow 
nonthermal linewidths < 0.1 kms -1 are among the most quiescent clumps ever 
found. The CCS clumps tend to exist around the higher density regions traced 
by NH-3 emission or submillimeter continuum sources, and the distribution is 
not spherically symmetric. Variation of the CCS abundance was suggested as 
an indicator of the evolutionary status of star formation. However, we can only 
find a weak correlation between iV(CCS) and n H2tVir . The velocity distributions 
of CCS clouds reveal that a systematic velocity pattern generally exists. The 
most striking feature in our data is a ring structure in the position-velocity 
diagram of L1544 with an well-resolved inner hole of 0.04 pc x 0.13 kms" 1 and 
an outer boundary of 0.16 pc x 0.55 kms -1 . This position- velocity structure 
clearly indicates an edge-on disk or ring geometry, and it can be interpreted as 
a collapsing disk with an infall velocity ^0.1 kms -1 and a rotational velocity 
less than our velocity resolution. Nonthermal linewidth distribution is generally 
coherent in CCS clouds, which could be evidence for the termination of Larson's 
Law at small scales, ~ 0.1 pc. 

Subject headings: ISM: abundances, ISM: clouds, ISM: kinematics and dynamics, 
ISM: molecules, radio lines: ISM, stars: formation 
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1. Introduction 

The CCS molecule was first discovered toward TMC-1 by Suzuki et al. (1984) and was 
identified in the laboratory by Saito et al. (1987). It is a very interesting molecule because 
it is abundant in starless clouds (Suzuki et al. 1992, hereafter SYO) with an appropriate 
critical density (10 4 — 10 5 cm -3 ) for pre-protostellar cores (Wolkovitch et al. 1997, hereafter 
WLG). CCS is also an excellent probe of kinematic structure of dark clouds because of its 
lack of hyperfine splittings and its moderate optical depth. Moreover, the CCS structure is 
Hund's case b and as such has a large Zeeman splitting factor of 0.70 Hz/xG _1 at 33 GHz 
(Shinnaga & Yamamoto, 1999), which provides an opportunity to explore the magnetic 
field structure in the early stage of star formation. 

Several single-dish surveys of CCS have been done in the past 10 years. The survey 
done by Fuente et al. (1990) showed that N(CCS) increases with N(C 18 0) and N(HC 5 N). 
However, this correlation is not so obvious for later observations (§4.2). SYO performed 
the largest single-dish survey toward 49 dark cores and concluded that CCS is more 
abundant in cold quiescent dark clouds while less abundant in star-forming regions. They 
also constructed a pseudo-time-dependent chemical evolution model, suggesting that CCS 
is more abundant in the early stage of chemical evolution whereas NH3 is more abundant 
in the later stage. [CCS]/[NH 3 ] was proposed as an indicator of cloud evolution and star 
formation. Based on this suggestion, Scappini & Codella (1996) observed NH 3 and CCS 
toward a sequence of objects in different evolutionary stages - dark clouds, inactive Bok 
globules, active Bok globules, and massive protostars. However, because of the large scatter, 
their data did not determine a trend. The most recent survey done by Benson et al. (1998) 
only showed a slight tendency for CCS to have a higher column density in starless cores, 
but their result was not statistically significant. 

CCS maps of dark clouds have been made toward B335 and L1498. Velusamy et al. 
(1995) combined VLA and DSN 70 m single-dish maps of the 22 GHz CCS J N = 2 1 - 1 
transition toward B335 and showed that the distribution of CCS emission is non-spherical 
and clumpy with depletion in the central 10". They concluded that CCS is characteristic 
of early astrochemical time (t< 10 5 yr) and more prominent in the region just outside the 
infall envelope. WLG showed that the CCS emission of L1498 has two peaks separated by 
2' in the SE-NW direction. However, the NH 3 map ( |Benson fc Myers 1989| ) only shows a 
single component between two CCS components. Kuiper et al. (1996) interpreted this as 
the CCS originating in a shell outside of the ammonia core, and concluded that L1498 may 
be about to collapse rapidly to form a protostar. 

In this paper, we present a survey with largest number of CCS maps and spectra. 
These data with large field of view (~ 6') and adequate velocity and spatial resolutions 
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provide a better foundation to examine the existent arguments. In §2 we describe our 
observation and data reduction procedure. In §3 we present the analysis for obtaining 
physical parameters and discuss the results in §4. We summarize the conclusions in §5. 

2. Observations 

We used 1-cm receivers (Carlstrom 1996) on the Berkeley-Illinois-Maryland Association 
(BIMA, Welch et al 1996) array to observe the CCS Jjv = 3 2 — 2i line toward 11 dark clouds 
or star- forming regions. The sources we observed are listed in Table 1. The frequency of 
this CCS transition is 33.75138 GHz (WLG). 

Most of the observations were carried out between 1997 Aug 27 and Sep 3 with one 
additional track for extended Bl emission on 1998 Sep 7. During the observations in 
1997, BIMA was in the D configuration with 9 antennas. The range of the baselines of D 
configuration was 0.6-8.0 kilowavelength, and the synthesized beam was about 30"x30", 
slightly different for different sources. The primary beam was about 6' at 33.75 GHz. All 
correlators to the most distant antenna from the center of array were unstable during the 
observations in 1997. The corrupt baselines were carefully flagged out. The extra track 
on Bl was obtained in order to observe extended emission around the edge of the primary 
beam in the previous observation. In this track, a northern field, B1(N), and a southwestern 
field, Bl(SW), were observed for 15 min each sequentially. The distances from the phase 
centers of each field to the center of the previous observation were less than half of the 
primary beam. BIMA was in a special compact configuration for this observation. Baseline 
range was about the same as D configuration, but because the compact array had more 
shorter baselines, it produced a 43"x36" beam for Bl. 

In order to resolve the extremely narrow linewidth of the quiescent dark clouds, the 
correlators were set in the mode with the highest resolution available on BIMA, 0.054 
kms -1 . In this mode, we only had one spectral window in each sideband. The total 
bandwidth of each spectral window was 6.25 MHz (55.3 kms -1 ). 

The 1-cm receivers operated from 27 to 36 GHz with a single sideband response at 
the bottom and top of the tuning range and with double sideband response in the center. 
Because the lower sideband was 2-3 times more sensitive than the upper sideband, we set 
the observed line in the lower sideband and the upper sideband was out of the operating 
range. As a result, we obtained a single sideband system. The system temperatures ranged 
from 40K to 96K during the observations in 1997 and from 130K to 230K for the observation 
in 1998. The much higher system temperature in the second year was due to worse weather. 
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We chose the brightest quasars within 10° of the sources as the phase calibrators. The 
phase calibrators were observed for 5 minutes every 30 minutes on the sources. We also 
observed a planet for 5 minutes for each source in order to calibrate flux. Because the 
observed bandwidth was very narrow and the signal to noise ratios were not high for the 
short observation time for each source, passband calibration was unnecessary. 

All data reduction was carried out using the MIRIAD software (Sault, Teuben, & 
Wright 1995). Bad visibility data were flagged, and data were phase calibrated and flux 
calibrated before being transformed into images. For those sources with line detections, 
the continuum was subtracted in the u — v plane, although the continuum was within 
the noise level. Natural weighting was used for constructing the images, and the CLEAN 
algorithm was used to deconvolve the dirty beam from the image. The boundary of the 
regions for deconvolution were chosen approximately between -la and 1 a rms levels 
around the sources in dirty maps, which provided better images especially for sources with 
high-sidelobe dirty beams. The rms noise level of all sources in Table 1 were estimated 
from a region in the clean map which did not contain signal. Primary beam correction was 
applied to the three fields of Bl which were then combined linearly. Taper was applied on 
the combined map of Bl to avoid excessive noise amplification at the edge of the mosaiced 
field and to achieve approximately uniform noise across the image. 



3. Results and Analysis 

Out of the 11 objects, 7 are detected with strong CCS emission and extremely 
narrow linewidthes (0.1-0.8 kms" 1 ). Among the 4 objects without CCS detection (DR21, 
OMC-N4, S106, W3), all except OMC-N4 are active star-forming regions, which suggests 
that CCS could be destroyed in the later stage of star formation. OMC-N4 is a dense 
(~ 10 6 cm" 3 ) region at the position north of the OMC1 core without active star formation 
where N bearing molecules are especially abundant. The 7 clouds with CCS detection show 
tremendous clumpy structure. We identified the clumps in these CCS clouds and analyzed 
them separately. A clump is identified in channel maps if there is an isolated peak with 
intensity greater than 3 a in the image-velocity data cube across more than 3 channels. 
We treat B335 and L63 as one clump each. B335 has double peaks, but it was previously 
modeled as a collapsing cloud (Zhou et al. 1993). L63 has the lowest signal to noise ratio, 
although it may separate into two clumps at the 2 a level. 

The peaks of individual clumps are listed in Table 2, and channel maps (Fig. l-7(a)) 
show the clumps marked with a letter designation. The Vi sr and linewidth of each clump 
are fitted with a Gaussian at the peak position and are listed in Table 2. 68% of clumps 
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have total linewidths < 0.25 kms -1 . If we exclude Bl and B335, both containing point-like 
IRAS sources, only one clump out of 16 has a total linewidth > 0.25 kms -1 . Assuming the 
kinetic temperature T& = 10 K, 7 clumps with nonthermal linewidths < 0.1 kms^ 1 are 
among the most quiescent clumps that have been found. The integrated intensity maps are 
shown in Fig. l-7(b). Fig. 1(c) shows the central velocity map of Bl, and Fig. 2-7(c) show 
interesting position-velocity diagrams for the rest of clouds. The distribution of CCS is 
discussed in §4.4. Fig. 8 shows the line profile of each clump. 

The brightness temperature of the intensity peak of each clump in Table 2 is obtained 
from the Rayleigh- Jeans Law, 



T b = 13.6 x ^Jy/bearn) M im ^ (1) 

where 6 a and 9 b are the dimensions of the synthesized beams listed in Table 1. To calculate 
the optical depth and the total column density, we need the excitation temperature T ex of 
CCS. The multiple-line CCS survey of SYO determined T ex = 5K for Bl, B335, L63, and 
L1544, and T ex = 5.2 K for L1498. We adopt T ex = 5K for all our samples except L1498. 
Because of the high intensity of L1498, its excitation temperature must be greater than 
5.5 K unless the optical depth is negative. We adopt a value T ex = 6K from Benson et 
al. (1998). Therefore, the optical depth r in Column 5 of Table 2 was obtained from the 
radiative transfer equation 



T = -ln(i - J(T b )/J(T ex )), (2) 

where J(T) is the Plank function. The total column density of CCS in Column 6 of Table 
3 was calculated from 



N(CCS) = 2.37 x 10- 9 gggj ex P[ E f Q T ^ W (cm - 2) , ( 3 ) 
K ' (2 J +1) a b A K J 

where Q(T ex ) is the partition function, E\ is the energy of lower level, A is the Einstein 
A-coefncient for the transition, and W — J l v dv is the integrated line intensity. We adopted 
frequencies of transitions and the Einstein A-coeflicient from WLG : Q (5K) = 23.8017, Q 
(6K) = 30.5570, E t = 1.6061 K, and A = 1.5999 x lO^s -1 . 

The virial masses in Column 7 of Table 3 were determined by assuming a self-gravitating 
sphere with n oc r~ l and a gaussian velocity distribution for each clump. MacLaren et al. 
(1988) derived the virial masses of such clouds as 
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M mr = 190 RAVI ( M ©)> ( 4 ) 

where R is the cloud radius in pc and AVh 2 is the line-of-sight FWHM speed of H 2 
molecules in kms -1 . If the density law is between n oc r^ 2 and n oc constant, the virial 
mass would only change from 0.6 to 1.1 M vir as given in equation (4). The slope of the 
iV(CCS)-n^ 2i „j r correlation derived in §4.2 is independent of the assumed density law. The 
radii of the clumps R listed in Columns 4 and 5 of Table 3 are estimated from half of the 
geometric average of the FWHM of major and minor axes of the clumps. Assuming the 
nonthermal part of the CCS and H2 linewidth are the same, AVh 2 can be calculated from 
AVccs an d the difference between their thermal widths. 



hT, 

AVl = 8ln(2)^ + AV 2 T 
m H2 



AV* cs + 8ln(2)kT(— — ), (5) 

inH2 mccs 



where is the kinetic temperature. We adopt Tk—10 K for all clouds. Because the observed 
AVccs is only a few channels wide for some clumps, the correction for the convolution effect 
is necessary. The true AVccs is estimated from {AVccs, true) 2 — (AVccs,obs) 2 — 0.054 2 . If 
we assume all the clumps are in viral equilibrium, we can calculate the H 2 density from the 
virial mass and estimate the fractional abundance of CCS, X(CCS) = N C cs/(2R ■ n H2iVir ). 
The results are listed in Column 8 and 9 of Table 3 and range from 0.1-9 xlCT 10 , which 
brackets the value 4.9 x 10 -10 obtained from the LVG calculation for L1498 by WLG. Large 
uncertainties for M vir , H 2 density, and the fractional abundance come from the unknown 
dimension along the line of sight. If the estimated dimension is larger by a factor of 2, then 
the M V i r would be larger by the same factor, the H 2 density would be smaller by a factor of 
4, and the fractional abundance of CCS would be larger by a factor of 2. 



4. Discussion 
4.1. Individual Clouds 

4.1.1. Bl 



Bl is an extensive molecular cloud in the Perseus complex. A distance 350 pc is 
adopted for comparison with other results, although there is some possibility that Bl may 
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be located in the other molecular component at less than 200 pc along the line of sight 
toward Perseus flCernicharo et al. 1985| ). Physical properties of the envelope and main core 



of Bl have been determined by the multiple-line study of Bachiller et al. (1990). In their 
study, the main core traced by NH3 has a linewidth of 0.8 kms -1 and a density greater 
than 10 4 cm" 3 . There are two IRAS sources (IRAS 03301+3057 and IRAS 03304+3100) in 
the observed fields (Fig. 1). These two source may play a role in destroying the physical 
condition for CCS to exist. 

We identify eight clumps in the channel map (Fig. 1). Clump A-F correspond to 
the NH3 main core, and the CCS distribution of the main core is very similar to the NH3 
distribution observed by Bachiller et al. (1990). The total mass of the main core estimated 
from the sum of the virial masses of clump A-F is 13.6M , which is much less than that 
derived from NHa(~ 150M©) ( |Bachiller et al. 1990| ). This can be understood if NH3 



traces the densest region of the core and CCS traces the envelope of the core (see §4.4). 
A systematic velocity field with a range from 6.0 kms -1 to 7.1 kms -1 is found along an 
arc formed by clumps A, B and F. This appears to be a rotational motion about IRAS 
03301+3057. The arc also has larger linewidths than that of other clumps in Bl, and the 
linewidths of all clumps are smaller than the width of NH 3 lines. B1(N) and Bl(SW) both 
have almost coherent central velocity. The large velocity differences between them and the 
nearest ends of the main core indicate that these two clumps are kinematically distinct from 
the main core system. 



4.1.2. B133 

B133 is a Bok globule at a distance of less than 400 pc; it has an IRAS source, IRAS 
19035-0657, at ~2.5' east of the center. Four clumps can be identified in the peak intensity 
map (Fig. 2). Clump A is weaker and has a narrower linewidth than the other clumps. 
Clumps B, C, and D may spatially connect to each other, and clump A may be isolated 
from them. Clump B has double peaks, which may be an infall feature or just simply two 
velocity components along the line of sight. The position- velocity diagram along C-B2-D 
shows a slight velocity gradient, about 1 ms _1 arcsec _1 = 2.5xl0~ 3 ms _1 AU _1 . The 
FWHM also increases from C, B, and to the southern east edge of clump D. 
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4.1.3. B335 

B335 is a low-mass star-forming Bok globule at a distance of 250 pc. A cold, compact 
far- infrared source lies in its center ( [Keene et al. 1983| ) , and it also shows dynamical features 
of a collapsing cloud with a bipolar outflow ( Frerking fc L anger T982D and infall (|Zh ou et 



|al. 1993|) . The observed asymmetric line profiles have been modeled as inside-out collapse 
( [Zhou et al. T993] , |Choi et al. 1995ft . Velusamy et al. (1995) combined VLA and DSN 70 m 



single-dish observations of the 22 GHz CCS Jn = 2i — lo transition and show that the 
distribution of CCS emission is non-spherical and clumpy with depletion in the central 10". 
They conclude that CCS is characteristic of early astrochemical time (t< 10 5 yr) and more 
prominent in the region just outside the infall envelope. 

The source shape in our integrated intensity map of B335 is similar to that of the DSN 
70 m map of Velusamy et al. The central position of the collapsing region (|Zhou et al.| 



1993| ) is offset toward the east of the integrated intensity peak by half of our synthesized 
beam. We obtain an average X(CCS)~ 1.6 x 10~ 10 , which is consistent with the values 
from the LVG calculation by Velusamy et al. : X(CCS)~ 3 x 10~ 9 in the outer envelope 
and X(CCS)~ 5 x 10~ u in the central depletion region. The dimension of the region with 
blue asymmetric line profiles can be estimated from the position-velocity diagrams along 
the major and minor axes; it is ~ 100"x30", corresponding to 0.12 pc x 0.04 pc, which is 
larger than the CCS depletion region. 



4.1.4. L63 

L63 is an isolated globule at a distance of 160 pc. A GO type star is located 
~2.5' southwest of the center. If this star is associated with the cloud, it may explain why 
L63 is the weakest CCS source in our sample. Two condensations at different velocity can 
be separated in the position- velocity diagram. The western one is centered at 5.8 kms -1 , 
and the eastern one at 5.9 kins" 1 . There is a velocity gradient along the east-west direction 
~ 1.4 ms _1 arcsec _1 = 8.7 xl0~ 3 ms 4 All -1 . The linewidth of the peak we obtained, 
0.167 kms -1 , is narrower than that obtained by SYO (0.8 kms" 1 ) and Benson et al. 1998 
(0.36 kms" 1 ). 



4.1.5. L183 



L183, or L134N, is a part of the L134-L183-L1778 system of interstellar clouds at high 
galactic altitude, b~ 36°. No far-infrared source has been detected toward L183 ( [Sargent et 
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al. 1983). Clark & Johnson (1981) observed OH 6 cm and H 2 CO 2 cm lines for this system, 
and suggested the retrograde rotation of these cores can be explained by magnetic braking. 

Three clumps are identified. Each clump is in a different velocity range. The flux of 
clump A could be stronger than the value in Table 2, because the intensity peak of clump 
A lies on the edge of the primary beam. Clump B has a different shape in channel maps, 
and it may divide into more clumps if the signal to noise ratio of the map were higher. 
The velocity gradient along A, B is obvious in the position- velocity diagram, with value 1.2 
ms 4 arcsec~ 1= 7.5 xlO -3 ms _1 AU _1 . 



4.1.6. L1498 

L1498 is an extremely quiescent pre-protostellar core in the Taurus complex. It has 
been studied extensively with multiple transitions of CCS at 22 GHz, 45 GHz, and 94 GHz 
observed by WLG. Our L1498 map clearly shows two concentrations in the position-velocity 
diagram along the major axis, which is consistent with the map obtained by WLG. The 
peanut-like distribution in the p-v diagram shows the two concentrations are at about the 
same velocity of 7.87 kms -1 and have a similar linewidth of 0.18 kms -1 . The virial mass of 
each clump is ~ 1.4 M , and the H 2 density derived from the virial mass is ~ 3xl0 5 cm -3 . 
The fractional abundances of CCS are 1.2 and 1.5xlCT 10 for the NW and SE clumps, and 
both values are less than that obtained from the LVG calculation by WLG , 4.9 xlCT 10 . 

4.1.7. L1544 

L1544 is a starless core in the Taurus complex at 140 pc. It is not associated with any 
point-like source. IRAS 05013+2505, peaked at 100 /im, is nearby, but is extended over 5'. 
Tafalla et al. (1998) did a multiline study of L1544 and found that a large area in L1544 has 
a double peak spectrum. They explained this feature as an infall with high self-absorption 
in the center. They modeled the velocity structure with a simple two-layer collapsing model 
and concluded that the data are inconsistent with the inside-out collapse model. They also 
ruled out ambipolar diffusion driven collapse because the speed of ionized species is similar 
to that of neutrals. 

Six clumps are identified in LI 544, which is elongated along the northwest to southeast 
direction. The position-velocity digram along the major axis of this elongated cloud shows 
a ring structure with an inner hole of 0.04 pc x 0.13 kms -1 and an outer boundary of 
0.16 pc x 0.55 kms" 1 . This ring structure can be interpreted as a collapsing disk or ring 
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almost edge-on with infall velocity ^0.1 kms -1 . The rotation velocity can be estimated 
from the velocity difference between the upper and bottom edge of the ring, which is only 
~ 0.05 kms -1 , comparable to our spectral resolution. Because the self-absorption toward 
the center has been observed with CS, the central hole of the ring should not indicate a 
physical hole in the center. Instead it may be a further evidence that CCS only traces the 
outer envelope of dense cores. 

4.2. Correlation of iV(CCS) with Density 

It has been widely suggested that CCS is more abundant in the early stage of 
protostellar cores. However, observationally it is difficult to determine a clear inverse 
correlation of iV(CCS) with H2 density as an indicator of the cloud age. Except for the 
work done by Fuente et al. (1990), observational data show that the correlation is fairly 
weak (SYO, WLG, Benson et al. 1998). Fuente et al. (1990) show a positive correlation 
between iV(CCS) and N(C 18 0) for seven dark clouds, which seems inconsistent with the 
prediction of chemical models unless the seven clouds are all in a very early stage. 

Since CCS distributions are clumpy and may be shell-like around the higher density 
core whereas CO apparently exists in a larger area, the depth of the CCS column should be 
a few factors smaller than that of the CO column. Therefore, the average densities obtained 
from CO may not be able to represent the density of CCS clumps. In addition, based on 
the complex distributions and active chemical reactions ( formation in the outer lower 
density region and destruction in the central higher density region), it seems unsuitable to 
use a single average density to represent the chemical behavior of CCS in the whole cloud. 
Therefore, it is necessary to analysis the iV(CCS)-n(H 2 ) correlation at the clump level. 

In order to have a systematic analysis at the clump level, we adopt the virial density 
to represent the H 2 densities. If the clouds are not virialized, the H 2 densities will be 
incorrect, but the statistical conclusion would still be valid if the virial densities were 
scaled from the actual densities with similar factors. We obtain a very weak correlation, 
iV(CCS) - io 13 - 7±LO n H ° 2 - 25±0 - 18 cm" 2 (Fig. 9). This relation brackets the weak correlation, 
iV(CCS) ~ 10 14 - 1 n]f 2 22 cm -2 obtained by WLG for one of the three velocity components of 
TMC-1D. In their study, the correlations for the other two components of TMC-1D are too 
weak to be distinguished from a flat distribution. Therefore, we conclude that even though 
CCS is more abundant in the early stage of a contracting cloud, the abundance of CCS is 
not likely to trace the age of the contracting cores very sensitively. Moreover, an almost flat 
distribution could imply that the time for creating and destroying CCS may be very short 
compared to the lifetime of CCS. 
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4.3. Cloud Equilibrium 

We compare the average H 2 column density obtained from C 18 or NH 3 observations 
(Myers et al. 1983, Benson & Myers 1989) to the average virial H 2 column density of each 
cloud for our CCS data. The average virial column density of each cloud is calculated by 
averaging that of all clumps in the cloud. The comparison results are listed in Table 4. 
Except for cloud Bl, the observed column density is much smaller than the CCS virial 
column density, which suggests that the CCS clumps are not virialized. Since the real 
densities of the CCS clumps should be smaller than the virial densities, X(CCS) estimated 
in Table 3 can only provide lower limits. 

4.4. Distribution of CCS 

The channel maps of our sampled clouds (Fig. l-7(a)) show the remarkable clumpy 
structure of the CCS distribution. As suggested by Velusamy et al. (1995) from the 
observation of B335, the clumpy structure may originate from the destruction of CCS in 
the densest region of the pre-protostellar cores. Therefore, combined with our conclusion 
in §4.3, we suggest that the clumps traced by CCS may not be self-contracting cores for 
protostars, but only higher density condensations in the envelope of a collapsing core. 

The above idea can be further supported by comparing the relative positions of CCS 
and higher density regions, traced by NH 3 emission or submillimeter continuum sources. 
We marked NH3 peaks with X symbols and marked submillimeter continuum sources with 
ellipses in the integrated intensity maps (Fig. l-7(b)). The comparison results are listed 
in Table 4. We found that except for cloud Bl, CCS generally exists around the higher 
density regions. This can be seen obviously for B133, L63, and L183, where the CCS 
emission exists asymmetrically around the submillimeter continuum sources. For L1498, 
the NH 3 peak appears between two CCS clumps, and for L1544, the NH 3 peak might be 
in the center of the ring structure discussed in §4.1.7. The NH 3 peak in B335 is close to 
the integrated intensity peak, and close to the place with central depletion reported by 
Velusamy et al. (1995). Although for the Bl main core (clump A-F) the CCS distribution 
is similar to the NH 3 distribution observed by Bachiller et al. (1990), the CCS distribution 
appears to be more clumpy than the NH 3 distribution. Also, the NH 3 peak is at the gap 
between clump B and E. Therefore, the CCS distributions for all clouds are consistent with 
the idea that CCS traces the envelope of pre-protostellar cores rather than being individual 
self-gravitation cores themselves. 
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4.5. Linewidth- Radius Relations 

Almost 2 decades after Larson (1981) discovered a clear correlation between the 
linewidths and the cloud sizes Av oc i? a38 , the physical origin of this Larson's Law 
remains an important unsolved problem. It has been widely discussed that turbulence 
or magnetohydrodynamic turbulence may provide the large nonthermal component of 
linewidths and the scaling relation through energy dissipation, though a complete theory 
has not yet been established. 

A question has been asked: does the relation terminate at small scale? Goodman et al. 
(1998) proposed a scenario that a transition to coherent linewidths takes place at the length 
scale ~ 0.1 pc based on their observations of NH 3 lines with spatial resolution 60"-80". The 
nonthermal components of linewidths in coherent cores are less than the thermal widths of 
H2, but strictly non-zero. They suggest that the coherent length scale may come from the 
dissipation threshold of MHD waves. 

The CCS clouds in our sample with physical sizes comparable to the coherent length 
scale can be a test ground of the existence of the coherent cores. The velocity dispersion of 
each pixel can be easily calculated from the second-order moment analysis. A cut-off level 
of 2 a on the spectrum must be set when calculating the velocity dispersion in order to 
reduce the confusion from noise. However, this high cut-off level can significantly reduce 
the velocity dispersion if the signal-to-noise ratio is low. Assuming an ideal Gaussian 
profile, we calculated how the correction coefficient varies with the signal-to-noise ratio. 
Averaging the corrected linewidth maps, we found that the linewidth-radius relation for 
each cloud is very consistent with a flat distribution. The result is not sensitive to the 
center position we chose. An example of L1498 is shown in Fig. 10. Since the linewidths are 
almost homogeneous everywhere for each cloud, the widths from the Gaussian fitting listed 
in Table 2 can represent the typical linewidth in the cloud. The total linewidth of CCS is 
almost contributed by the nonthermal width, since CCS is a heavy molecule. The thermal 
linewidth of H 2 at T^ =10 K is about 0.5 kms" 1 . Therefore, the nonthermal linewidths of 
our sampled clouds are comparable or less than the thermal linewidth of H2. 

The observed facts seem to support the picture of coherent cores. However, as we 
discussed in §4.4, we have to be cautious that CCS may only trace the envelope rather than 
the center of cores. Therefore, the linewidth at small radius on the sky does not necessarily 
represent the linewidth at small distance from the center of the core. Nevertheless, we can 
say at least the linewidths are very homogeneous in the narrow density regions traced by 
CCS. 
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5. Summary 

CCS is an abundant molecule in quiescent dark clouds and is not detectable in active 
star-forming regions. Seven out of 25 clumps have extremely small nonthermal linewidths, 
< 0.1 kms -1 . The distributions of CCS are very clumpy, and the CCS clumps are not 
necessarily symmetric around the higher density regions where protostars may be born. 
The CCS clumps themselves are not likely to collapse into protostars, but rather exist 
transiently in the regions with appropriate densities as the cloud evolves. Therefore, we 
conclude that CCS traces pre-protostellar cores in the very early stage of star formation or 
traces the envelope of pre-protostellar cores in the later stage. Homogeneous linewidths in 
CCS clouds may provide evidence that the nonthermal linewidths transition to coherence 
at the length scale ~ 0.1 pc, which is proposed by Goodman et al. (1998). Because CCS is 
a Zeeman molecule and also an excellent kinematic probe, follow-up polarimetry and high 
velocity resolution studies could provide a good opportunity to understand the physical 
properties of star-forming regions in the early stage. L1544, with a very systematic velocity 
field and strong CCS emission, is an excellent target to model the physical condition of a 
collapsing core. 

Research at the Laboratory for Astronomical Imaging was supported by NSF grant 
AST 96-13999 and by the University of Illinois. 
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Fig. I. — (a) CCS channel maps of Bl obtained from 3 mosaic fields. The rms noise 
level varies across the maps and ranges from 0.1 Jybeam -1 at the center of the central 
field (R.A.(J2000)=03:33:17.83, Dec.(J2000)=31:07:29.80) to the maximum value of 0.19 
Jybeam -1 beyond 3' from the center. Contours are drawn at -0.5, -0.3, 0.3, 0.5, 0.8, 1.2 
Jybeam -1 . The black square and star indicate the positions of IRAS 03301+3057 and 
03304+3100 (=LkHo;327). (b) The integrated intensity map of Bl. Contours are drawn at 
-0.35, -0.25, -0.15, 0.15, 0.25, 0.75 Jybeam -1 kms -1 . The X symbol marks the NH 3 peak 
(Bachiller et al. 1990). (c) The intensity weighted velocity of Bl. Contours are at 6.2, 6.5, 
6.8, 7.1 kms -1 . 

Fig. 2. — (a) CCS channel maps of B133. The 1-a rms noise level is 0.33 Jybeam -1 , and 
contours are drawn at -3, 3, 5, 7, 9 a levels, (b) The integrated intensity map of B133. 
The 1-a rms noise level is 0.1 2 Jy beam -1 kms -1 , and contours are drawn at -3, 3, 6 a 
levels. The X symbol marks the NH 3 peak (Benson & Myers 1989), and the circle marks 
the submillimeter continuum source (Ward-Thompson et al. 1994). (c) The position- velocity 
diagram along the cut shown in (b). Contours are at 15%-90% of the peak, and the contour 
interval is 15%. 

Fig. 3. — (a) CCS channel maps of B335. The 1-a rms noise level is 0.22 Jybeam -1 , and 
contours are drawn at -3, 3, 5, 7, 9 a levels. The dot marks the position of the spectrum 
shown in Fig 8. (b) The integrated intensity map of B335. The 1-a rms noise level is 0.10 
Jybeam -1 kms -1 , and contours are drawn at -3, 3, 5, 7 a levels. The X symbol mark the 
NH 3 peak (Benson & Myers 1989), and the plus symbol marks the center of CCS depletion 
(Velusamy et al. 1995). (c) The position- velocity diagram along the cut shown in (b). 
Contours are at 30%-90% of the peak, and the contour interval is 15%. 

Fig. 4. — (a) CCS channel maps of L63. The 1-a rms noise level is 0.22 Jybeam -1 , and 
contours are drawn at -2, 2, 3, 5 a levels. The dot marks the position of the spectrum 
shown in Fig 8. (b) The integrated intensity map of L63. The 1-a rms noise level is 0.035 
Jy beam -1 kms -1 , and the contour is drawn at 3 a level. The X symbol marks the NH 3 
peak (Benson & Myers 1989), and the ellipse marks the submillimeter continuum source 
(Ward-Thompson et al . 1994). (c) The position- velocity diagram along the cut shown in 
(b). Contours are at 30%-90% of the peak, and the contour interval is 15%. 

Fig. 5.— CCS channel maps of L183. The 1-a rms noise level is 0.37 Jybeam -1 , and 
contours are drawn at -3, 3, 5 a levels. The plus symbol marks the pointing center, and 
clump A is on the boundary of the primary beam, (b) The integrated intensity map of 
L183. The 1-a rms noise level is 0.035 Jy beam -1 kms -1 , and contours are drawn at -6, -3, 
3, 6, 9, 12 a levels. The X symbol marks the NH 3 peak (Benson & Myers 1989), and the 
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ellipse marks the submillimeter continuum source (Ward-Thompson et al. 1994). (c) The 
position- velocity diagram along the cut shown in (b). Contours are at 30%-90% of the peak, 
and the contour interval is 15%. 

Fig. 6.— (a) CCS channel maps of L1498. The l-a rms noise level is 0.15 Jybeam -1 , 
and contours are drawn at -5, -3, 3, 5, 8, 11, 14, 17, 20, 23, 26 a levels. The dots mark 
the positions of the spectrum shown in Fig 8. (b) The integrated intensity map of L1498. 
The l-a rms noise level is 0.050 Jybeam -1 kms -1 , and contours are drawn at -6, -3, 3, 6, 9, 
12, 15, 18 a levels. The X symbol marks the NH 3 peak (Benson & Myers 1989). (c) The 
position- velocity diagram along the cut shown in (b). Base contour is at 3-a rms noise level, 
0.45 Jybeam -1 , and the contour interval is also 0.45 Jybeam -1 . 

Fig. 7. — (a) CCS channel maps of L1544. The l-a rms noise level is 0.13 Jybeam -1 , and 
contours are drawn at -6, -3, 3, 6, 9, 12, 15, 18 a levels, (b) The integrated intensity map 
of L1544. The l-a rms noise level is 0.050 Jybeam -1 kms -1 , and contours are drawn at -3, 
3, 6, 9, 12 a levels. The X symbol marks the NH 3 peak (Benson & Myers 1989). (c) The 
position- velocity diagram along the cut shown in (b). Base contour is at 2-a rms noise level, 
0.26 Jybeam -1 , and the contour interval is also 0.26 Jybeam -1 . Note that the central two 
contours are at 0.26 and 0.52 Jybeam -1 , which indicate a hole in the center. 

Fig. 8. — Line profiles of CCS clumps. 

Fig. 9. — CCS column density vs. virial H 2 density for all clumps (see §4.2). Note that 
BID is an unresolved clump and L183A is on the edge of the primary beam. 

Fig. 10. — The linewidth-radius relation of L1498. The central position is chosen at the 
intensity peak. The rise feature at R>0.1 pc can be washed out by choosing different central 
positions, however, the flatness within R=0.1 pc remains unchanged. 
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Table 1. Observation parameters 
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Source 


R.A.(J2000) 
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B335 
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32 


24 
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Y 
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0.21 


38 


30 


71 
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L183 
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22 
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37:22:54.00 
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0.23 
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0.19 


39 
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Table 2. Line parameters 



Source 


Distancc a 
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-* peak 
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(K) 
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(km/s) 


B1A 
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1.36 


1.61 


0.48 


6.11±0.01 


0.42±0.03 


BIB 




0.73 


0.87 


0.23 


6.78±0.03 


0.78±0.08 


B1C 




0.81 


0.96 


0.26 


6.48±0.02 


0.24±0.04 


BID 




0.70 


0.83 


0.22 


6.50±0.01 


0.13±0.03 
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1.21 


0.34 
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0.61±0.06 
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1.23 


1.46 


0.42 


6.97±0.01 


0.38±0.03 


Bl(SW) 




2.12 


2.51 


0.90 


6.81±0.01 


0.34±0.02 


B133A 


400 


1.87 


1.90 


0.61 


11.96±0.01 


0.13±0.03 


B133B1 




3.03 


3.14 


1.4 


12.10±0.01 


0.15±0.03 


B133B2 




3.16 


3.28 


1.5 


12.27±0.01 


0.09±0.02 
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2.43 


2.52 
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2.54 


2.63 
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L183A 
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a Reference: Hilton & Lahulla (1995). 

h T ex = 6K for L1498; T ex = 5K for the rest of clouds. 
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Table 3. Derived quantities 
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Table 4. Comparison between CCS clouds. 
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